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基于 策略 改进 遗传 算法 的 双星 星 族 光谱 拟 合 
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671000; 3. 广州 大 学 天 体 物 理 中 心 ， 广 东 广 州 510006) 


摘 要 :， 双 星星 族 光谱 拟 合 的 实质 是 在 双星 星 族 理论 光谱 库 对 应 的 参数 空间 进行 搜索 ， 
通过 最 小 二 乘法 找到 与 被 拟 合 观测 光谱 最 相近 的 理论 光谱 。 快速 、 准 确 的 拟 合 是 有 效 利 用 星 
族 合成 方法 处 理 海量 星系 光谱 的 难点 和 关键 . 为 了 提高 双星 星 族 光谱 拟 合 的 速度 , 针对 双星 
星 族 光谱 拟 合计 算 量 巨 大 、 拟 合 效 率 不 高 等 问题 ,本 文 在 双星 星 族 光 谱 拟 合 的 优化 模型 基础 
上 ,采用 策略 改进 的 遗传 算法 对 模型 进行 求解 ,并 与 BS2fit 算法 和 传统 遗传 算法 进行 比较 。 
实验 表明 ,通过 改进 策略 的 遗传 算法 对 双星 星 族 光谱 拟 合 的 速度 平均 可 提高 43. 5h, 一定 程 
度 上 推动 了 演化 蜂 族 合成 方法 在 天 文 研究 中 的 应 用 。 
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双星 星 族 是 包含 一 些 共同 演化 的 双星 和 单 星 , 具有 相同 初始 金属 丰 度 和 演化 年 龄 的 恒星 
组 合 , 它 是 对 单 星星 族 的 扩展 和 丰富 。 基 于 双星 星 族 模 型 得 到 的 星团 和 星系 参数 结果 通常 更 
合理 "， 因 此 ， 基 于 双星 星 族 的 研究 得 到 越 来 越 多 的 关注 。 

光谱 拟 合 是 获得 星系 物理 信息 的 重要 方法 ”。 随 着 各 种 巡天 项 目的 开展 和 天 文 信息 技术 
的 发 展 ， 人 们 获得 了 越 来 越 多 的 高 质量 光谱 数据 ， 为 天 体 光 谱 拟 合 创造 了 条 件 ， 有 效 地 推动 
了 光谱 拟 合 的 发 展 。 光 谱 拟 合 近年 来 取得 了 一 系列 成 果 ， 如 文献 [2 成功 地 用 贝 叶 斯 方法 对 
时 星星 族 光谱 拟 合 估计 出 星系 参数 信息 。 文 献 [3] 基 于 贝 叶 斯 证 据 用 贝 叶 斯 模型 选择 的 方法 
分 析 了 双星 星 族 和 单 星星 族 模型 ,发现 双星 星 族 模型 比 其 它 没 有 考虑 双星 相互 作用 的 模型 的 
光谱 拟 合 效果 好 很 多 ,文献 [4] 和 [5] 对 早期 星系 紫外 到 可 见 波段 的 光谱 拟 合 发 现 双星 星 族 模 
型 对 球状 星团 蓝 离散 星 观测 结果 有 很 好 的 解释 。 文 献 [6]Han, Z. (2007) 首次 通过 双星 星 族 光 
谱 拟 合 发 现 它 能 很 好 地 解释 早 型 星系 的 紫外 反 转 现象 。 

由 于 双星 星 族 演 化 比 传统 单 星星 族 更 复杂 , 考虑 的 物理 输入 也 更 多 , 其 对 应 的 理论 光谱 
JERK, 一 定 拟 合 精度 要 求 下 的 光谱 采样 率 也 更 高 , 所 以 双星 星 族 光 谱 拟 合 需 要 进行 海量 的 
计算 。 能 否 快 速 高 效 地 进行 双星 星 族 光谱 拟 合 直接 影响 着 星 族 演化 合成 方法 的 应 用 。 因 此 ， 
探索 基于 双星 星 族 模型 的 快速 有 效 的 光谱 拟 合 方法 显得 十 分 重要 。 

为 了 解决 双星 星 族 光谱 拟 合 速 度 不 能 满足 星系 演化 等 研究 需求 的 问题 , 加 快 星 族 演化 合 
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成 方法 在 天 文 


化 为 一 类 最 优化 问题 ， 


1 双星 星 族 模型 介绍 


ChinaXiv 合 作 期 刊 


研究 中 的 推广 和 运用 , 本文 基于 管理 策略 进行 了 一 个 尝试 。 我们 将 光谱 拟 合 转 


以 最 小 化 理论 光谱 和 观测 光谱 之 间 的 均 方 误差 为 优化 目标 , 建立 最 优 
化 模型 ， 提 出 改进 策略 的 遗传 算法 对 模型 求解 ， 并 与 BS2fit 算法 "和 传统 遗传 算法 对 比 以 
验证 算法 的 有 效 性 。 


双星 星 族 模型 的 一 个 主要 目标 是 构建 高 精度 且 完 备 的 双星 星 族 理论 光谱 库 。 双 星星 族 的 


理论 光谱 构建 包括 两 个 部 分 :模型 输入 参数 的 假设 和 理论 光谱 的 计算 。 
如 文献 1] 所 述 , 双星 星 族 模型 采用 Hurley 


量 比 q 服 从 0 
示 星 系 的 年 龄 

利用 演化 
同时 刻 形 成 的 


的 快速 恒星 演化 程序 代替 恒星 演化 轨迹 , OK 


用 BaseL3. 1 恒星 光谱 库 和 Chapbrier 初始 质量 函数 ”"， 并 假设 每 个 双星 的 主星 和 次 星 质 
-1 的 均匀 分 布 “。 恒星 形成 率 假设 为 下 (t) = Y0) . era 的 形式 ,上 和 1 分 别 表 


和 恒星 形成 的 持续 时 标 。 


星 族 合成 方法 , 在 初始 质量 函数 和 恒星 形成 率 已 知 的 情况 下 可 以 计算 星系 中 不 
恒星 数 。 然 后 ， 借 助 Hurley 的 快速 恒星 演化 程序 来 进行 恒星 演化 计算 。 最 
通过 BaSeL3. 1 恒星 样本 库 将 恒星 的 演化 参数 转换 为 光谱 能 量 分 布 。 在 完成 恒星 演化 计算 


E] 
E] 


将 所 包含 的 单 颗 恒 星 的 光谱 进行 合成 ， 就 可 以 计算 年 龄 为 :的 星 族 在 波长 4 的 光谱 流量 有 0。 


每 个 星 族 光 谱 流 量 由 三 个 参数 决定 :金属 丰 度 Z, 年 老 星 族 年 龄 tt ， 和 年 轻 星 族 年 龄 tz。 考 虑 


其 中 ， fro 为 


FRB ICA, Ja iY 


三 个 参 数 (Z, ti, tz) 决定 星 
为 和 Go, 一 个 高 芭 


fh = [fao(z, ti t2) WO G(v,, 


能 量 比 ，s 为 卷 积 算 子 。 


为 了 简化 双星 星 
然而 这 并 不 影响 实验 目的 , 同样 可 
因为 我 们 所 采取 的 对 


恒星 速度 弥散 0, 的 影响 , 尘埃 消光 (采用 Calzetti 等 人 在 2000 年 提出 的 消光 理论 ”), 红 移 
银河 系 消光 等 因素 ， 双 星星 族 合成 中 在 波长 处 观测 的 光谱 流量 可 表达 为 : 


0.) p(t, )r (Ay) #(1) 


WRIGHT, G(v,, os) 是 平均 值 和 标准 偏差 分 别 
分 布 ,p (TX) 是 通过 尘埃 消光 为 的 尘埃 的 能 量 百 分 比 ,r (AQ) 是 通过 银河 


族 的 理论 光谱 库 , 本 文 将 b,,o,， Tr 4A4 的 取 值 均 设 为 0, 仅 考 虑 Z, ty, tzo 
以 对 比分 析出 所 提出 的 该 进 遗 传 算法 在 光谱 拟 合 中 的 表现 ， 


} 比 实验 也 是 基于 同样 的 简化 


2 基于 策略 改进 遗传 算法 的 双星 星 族 光谱 拟 合 


2. 1 双星 星 族 光 谱 拟 合 模型 


光谱 拟 合 


的 目的 是 为 了 获得 星系 星团 


的 双星 星 族 理论 光谱 库 。 


的 物理 参数 信息 。 如 上 一 节 所 述 , 在 双星 星 族 模型 


质量 可 由 Thomas 等 (2005) 9 的 拟 合 公式 计算 : 


式 中 的 M, 为 星 


型 
中 光谱 是 由 z, ty, tp, Vo Os Tn AX 和 星 族 恒星 质量 8 个 参数 共同 决定 的 。 其 中 ， 星 系 的 恒星 


log(M./M) = 0.63 + 4.5210g (00/ (kms*)) (2) 
系 的 恒星 质量 , oo 是 恒星 的 速度 弥散 。 其 它 参数 可 通过 光谱 拟 合 ， 即 将 星 族 的 


理论 光谱 与 观测 光谱 进行 比 对 从 而 给 出 星 族 组 成 


到 最 接近 观察 光谱 的 到 
小 化 理论 光谱 与 观测 六 


的 主要 参数 。 


的 定量 信息 。 具 体 任务 是 从 理论 光谱 库 中 找 


论 光 谱 。 在 光谱 拟 合 中 ， 通 过 将 观测 光谱 与 理论 光谱 进行 比 对， 并 最 
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6 谱 数 据 之 间 的 差异 以 确定 拟 合 最 优 的 理论 光谱 , 进而 推断 出 恒星 组 成 
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双星 星 族 理论 光谱 库 有 稠密 的 参数 网 格 : 恒星 金属 丰 度 (Z) 取 值 为 0. 00030. 03, 4 Z< 
0. 001 时 其 参数 间隔 为 0. 0001, 金属 丰 度 较 高 时 (Z > 0.001) 参数 间 隔 设 为 0. 001。 星 族群 
的 年 老 星 族 年 龄 (t1) 在 0 ”15Gyr 范围 内 , 其 间隔 设 为 0. 16yr ;年 轻 星 族 年 龄 tz) 从 t1 递 
减 到 0Gyr， 其 间隔 设 为 0. 1Gyr。 人 恒星 速度 弥散 6 ,参数 以 步 长 为 10kns 从 0 到 350km 
s 改变 ;尘埃 消光 (t) 取 值 是 时 间隔 为 0.01 的 0 至 0.5 区 间 。 因 此 ， 双 星星 族 理论 光谱 
库 包含 巨大 参数 空间 ， 交 谱 拟 合 包含 天 然 的 海量 计算 ， 对 此 建立 拟 合 模型 如 下 。 

记 观 察 光 谱 为 Fa ， 双 星星 族 合成 理论 光谱 为 Fa， 则 光谱 拟 合 契合 度 的 可 由 

X? = Lifa — fina) oa ]#(3) 


刻画 。 其 中 是 波长 的 权重 参数 。 最 佳 拟 合 结果 对 应 于 min (2) 。 在 大 多 数 情况 下 , wi 使 
用 观测 误差 B2 表示 : 


w, = 1/E2# (4) 


co tH FY BAAR Ba Re RE TERE. EE TBST > HS SE J EE VO SG ERN De HS 
间 进 行 搜索 ， 找 到 与 被 拟 合 的 观察 光谱 X 最 小 值 押 对 应 的 理论 光谱 。 当 X 最 小时, 双星 星 族 
合成 光谱 拟 合 模型 满足 物理 约束 : 0<w;<1， 其 本 质 是 一 个 最 优化 问题 ， 可 形式 化 为 : 


min x? 
s.t. wi>0 


Dwi=1 


其 中 ，(3) 式 为 目标 函数 。 
2.2 双星 星 族 光 谱 拟 合算 法 


遗传 算法 "…” (GA, Genetic Algorithm) 是 模拟 自然 界 生 物 遗 传 进化 行为 ， 提 出 选择 、 

交叉 和 变异 算 子 ,通过 对 求解 问题 空间 的 多 点 并 行 迭 代 搜 索 能 快速 有 效 的 收敛 至 全 局 最 优 解 
的 一 类 仿生 启发 式 算法 。 遗传 算法 是 一 种 基于 种 群 的 搜索 算法 , 其 个 体 代 表 了 可 能 解 的 一 个 
样本 。 随 着 搜索 的 迭代 进行 , 个 体 通过 分 享 和 交换 他 们 之 间 的 信息 来 产生 更 好 的 后 代 。 传统 
的 遗传 算法 (CGA, Conventional Genetic Algorithm) 包含 选择 、 交 叉 和 突变 等 至 关 重 要 
的 操作 。 他 们 可 以 选择 获胜 的 个 体 ， 维 持 种 群 的 多 样 性 ， 并 在 种 群 内 的 个 体 之 间 交 换 信息 。 

然而 ， 传 统 改 进 策略 遗传 算法 在 局 部 搜索 能 力 上 表现 较 差 , 且 对 初 值 非常 敏感 。 对 GA 
算法 运行 机 制 分 析 后 ， 受 粒子 群 算法 惯性 自 适 应 的 启发 ， 同 时 为 克服 该 不 足 ， 本 文 在 传统 遗 
传 算法 改进 策略 的 基础 上 设计 了 一 种 自 适 应 遗传 算法 (SGA, Self-adaptive Genetic 
Algorithm) ,来 求解 双星 星 族 光谱 拟 合 问 题 。 传 统 遗 传 算 法 (CGA) 采用 固定 交叉 概率 (PC, 
Probability of Cross) 和 固定 突变 概率 (PM, Probability of Mutation) 。 这 使 得 算法 
按照 既定 的 规则 进行 ,从 而 导致 遗传 变异 不 够 灵活 。 基 于 特定 的 进化 遗传 使 得 算法 的 性 能 如 
收敛 速度 和 全 局 优化 能 力 对 进化 参数 的 选择 很 敏感 。 如 果 个 体 的 适应 度 值 大 于 总 体 的 平均 适 
应 度 值 ， 则 分 配 的 交叉 概率 和 变异 概率 应 该 较 小 。 此 外 ， 对 适应 度 值 越 大 的 个 体 ， 表示 该 个 
体 越 接 近 目 标 解 , 因此 ,应 赋予 其 较 小 的 交 又 和 变异 概率 , 以 保护 其 信息 不 被 遗传 给 下 一 代 。 
与 之 相反 , 对 适应 值 较 小 的 个 体 则 应 伺机 引入 扰动 因子 增强 其 变异 能 力 , 个 体 应 被 赋予 较 大 
的 交叉 和 突变 概率 ， 以 增加 其 变异 概率 。 于 是 ， 适 应 度 值 越 小 的 个 体 变异 的 机 会 越 大 。 为 了 
减少 进化 参数 选择 的 依赖 性 ， 提 高 算法 的 收敛 速度 和 全 局 优化 能 力 ，SGA 改进 如 下 : 

1) KWN PC 更 新 : 
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me F - Lt)) 大 > fang 
- 


Co; fe < favg (5) 
0.3 
w= (Tea ) 
2) BEN PM 更 新 : 
_ (1-76a) (fa 一身 ) \. 
Pm = fag fmax— favg J Ím > favg 


6 
Mb; fn < fave (6) 


其 中 ，jiteration 为 当前 迭代 次 数 ， 所 ,为 种 群 的 平均 适应 值 ，f 为 选择 交叉 操作 的 个 
体 适 应 值 ， 所 为 选择 变异 操作 的 个 体 适应 值 ，fmax 为 最 大 适应 值 ，yos 为 随机 参数 ， 据 经 验 取 
Yoa = 0.5. 


D 干扰 因子 主要 目的 是 增强 算法 对 搜索 空间 的 探索 能 力 。 为 此 ， 引 入 了 如 下 扰动 
子 (DO, Disruption Operator) : 


Ba 


Dis;,; x 6(-2, 2) if Dis i, best 2 1 
Dy = 


1+ Disi, best X (2, 1) otherwise (7) 


其 中 : Sæ (-2, 2) 之 间 的 随机 数 ， Dis ij 是 第 i 个 染色 体 和 领域 第 j 个 染色 体 的 汉 明 距 
离 。Disi pest 是 当前 最 优 染 色 体 和 第 i 个 染色 体 的 汉 明 距离 。 


据 (6) 式 更 新 PM 


图 1 SGA 流程 图 


Fig. 1 SGA flow chart 
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SGA 算法 具体 流程 如 图 1 所 示 ， 主 要 包括 以 下 步 又， 

步骤 一 : 光谱 数据 预 处 理 , 包括 剔除 发 射线 、 归 一 化 和 红 移 处 理 ; 

步骤 二 :问题 编码 与 种 群 初始 化 ; 

步骤 三 : 按 目标 函数 计算 适应 值 并 保留 当前 最 优 种 群 。 决 定 种 群 的 大 小 包括 染色 体 长 
度 和 个 数 等 ， 依 策略 初始 化 种 群 ， 设 定 进化 终止 条 件 ; 

步骤 四 : 判断 是 否 达到 收敛 条 件 或 停止 条 件 ， 如 果 满 足 则 输出 搜索 结果 ， 否 则 进行 返 
代 进 化 ; 

步骤 五 : 选择 算 子 :根据 适应 度 值 随机 选择 一 组 较 好 的 个 体 ; 

WIR: 按 (5) 式 更 新 交叉 概率 ; 

步骤 七 : 交叉 算 子 :从 集合 中 随机 选择 两 染色 体 交 换 一 个 或 多 个 位 ， 根 据 交 叉 概 率 产 生 
两 个 新 的 个 体 ， 按 (7) 式 添加 扰动 ; 

步骤 八 : 按 (6) 式 更 新 变异 算 子 : 从 集合 中 随机 选择 一 个 染色 体 ， 并 根据 突变 概率 改 
变 一 个 或 多 个 位 以 产生 一 个 新 的 染色 体 ; 

步骤 九 : 返回 步骤 四 进行 终止 判断 ， 循 环 以 上 步骤 ， 直 到 终止 条 件 。 


3 实验 与 分 析 


3. 1 改进 策略 的 遗传 算法 实验 

为 了 验证 SGA 在 双星 星 族 光谱 拟 合 中 的 寻 优 效率 ， KAIR OOM LS 
算法 作为 对 比 。 实 验 代码 采用 Python3. 6 开发 ， 操 作 系统 为 ubuntul8. 04 64bit， 内 存 容 
8 GB, CPU 型 号 为 i7-4790, 主 频 3. 60GHz。 实 验 的 理论 光 训 sh Sic FSF Chapbrier il 如 
数 的 双星 星 族 理论 光谱 库 进 行 拟 合 。 观测 光谱 数据 选 自 99 个 本 地 星系 的 紫外 -可 见 波段 的 光 
谱 数 据 ， 可 以 在 互联 网 上 下 载 。 
这 些 数据 的 特点 是 覆盖 从 紫外 到 可 见 光波 段 的 波长 ， 范 围 很 广 。 此 外 ， 由 于 双星 对 紫 
外 光谱 有 明显 的 影响 , 这 些 数据 也 是 探索 双星 演化 对 光谱 拟 合影 响 的 理想 数据 。 每 条 观察 光 
谱 采 样 点 数 是 N=3871， 以 5500A 为 归 一 化 因子 ， 对 光谱 进行 预 处 理 。SGA 算 法 参数 设置 为 : 
交叉 概率 上 限 C2=1, 变异 概率 上 限 M2=1， 随 机 参数 yes = 0.5， 个 体 编码 长 度 n=20， 种 群 规模 
E Pop=100。 耗 时 指标 取 每 个 算法 重复 20 次 的 平均 值 。 
= 实验 首先 从 观测 光谱 数据 中 随机 选取 一 条 (NGC1399 星 系 ) 作 为 实验 数据 参考 文献 [12] 


© 进行 数据 预 处 理 ， 然 后 将 观察 光谱 进行 坏 点 、 发 射线 、 ee 
spec-0266-51630-0010 的 光谱 数据 ; 接着 ,借鉴 相关 研究 ,采用 x?/N, 作 为 评价 标准 ， 当 
77/7, 《1 时 ， 认 为 理论 光谱 与 观测 光谱 拟 合 良 好 ， BK or ia, 拟 合 结果 可 以 


接受 。 拟 合 迭 代 过 程 如 图 2 左 所 示 。 算 法 在 种 群 迭代 20 次 后 收 化 ， 对 应 的 适应 值 x? =1308, 
er Nt 理论 光谱 与 观测 光谱 拟 合 良 好 ， 如 图 2 右 所 示 。 其 中 灰色 曲线 
代表 实际 观察 的 光谱 ， 红 色 曲 线 代表 所 拟 合 得 的 最 优 理论 光谱 。 


!. http://cluster.shao.ac.cn/~zhongmu/ssp_chab_hrseds/ 


“. ftp://ftp.stsci.edu/pub/catalogs/nearby_gal/sed.html 
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— theory spectrum 
1335 A — observed spectrum 
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Iteration AIA 


图 2 SGA Waa ULE 


Fig. 2 SGA fit 


Co 


. 2 与 传统 遗传 算法 (CGA) 的 比较 
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0 1000 2000 3000 4000 5000 
lteration 
图 3 CGA 拟 合 迭 代 过 程 图 4 BS2fit 拟 合 遍历 过 程 
Fig. 3 CGA fit Fig.4 BS2fit fit 


CGA 算法 参数 : 交叉 概率 PC=0. 8, 变异 概率 PM=0. 1， 染 色 体 长 度 为 20，Pop=100。 对 测 
试 光谱 用 CGA 算法 在 双星 星 族 理论 光谱 库 上 进行 拟 合 结果 如 图 3 所 示 。CGA SRE, PEE 
代 16 次 后 陷入 局 部 最 优 , 迭代 (iteration)15 次 后 开始 收敛 , 寻 到 拟 合 度 最 优 的 理论 光谱 。 
SEK 20 次 平均 耗 时 3201. 3s 比 BS2fit 算法 略为 占 优 ， 速 度 提高 了 19. 2%。 与 CGA 算法 
相 比 ，SGA 在 迭代 20 次 后 收敛 ， 能 有 效 跳出 局 部 最 优 陷阱 ， 收 敛 速度 显著 改善 ， 实 验 20 次 
平均 耗 时 为 1304. 4s， 速 度 提高 了 19. 1%。 


3. 3 与 BS2fit 算法 的 比较 

BS2fit 算 法 采用 文献 [1] 所 述 的 BS2fit 算 法 python 实 现 , 拟 合 遍 历 过 程 如 图 4 所 示 , 横 
轴 表 示 遍 历 计算 的 次 数 ， 纵 轴 表 示 每 次 计算 得 到 的 卡 方 统计 量 的 数值 。 由 于 是 全 网 格 搜索 ， 
双星 星 族 模型 的 理论 库 履 盖 波 段 长 .演化 参数 网 格 稠密 (模型 光谱 库 履 盖 金 属 丰 度 Z 从 0. 0003 
到 0. 03, 精度 为 0. 0001;t1、t2 范 围 从 0 到 15Gyr, 精度 为 0. 1 Gyr) 等 导致 样本 参数 空间 巨大 。 
虽然 我 们 将 模型 参数 做 了 简化 处 理 , 但 是 要 完成 一 条 观测 光谱 的 拟 合 也 需要 逐一 进行 海量 对 
比 ， 能 寻 到 拟 合 度 最 优 的 理论 光谱 。 实 验 20 次 平均 耗 时 为 3954. 2s。 与 BS2fit 算 法 相 比 ，SGA 
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的 全 局 搜索 能 力 ， 对 理论 光谱 搜索 效率 更 高 ， 收 敛 速度 显著 改善 ， 20 次 平均 耗 


% 0 


如 上 ， 三 种 算法 均 能 找到 最 小 好 对 应 的 理论 光谱 ， 按 照 理论 光谱 的 参数 网 格 设置 ， 反 
演 得 到 表 1 中 给 出 的 拟 合 参数 。 


#1: NGC1399 星系 bsSP-fitted 拟 合 参 数 


Tab. 1 bsSP-fitted Parameters of NGC1399 Galaxies 


星系 参数 Z (金属 丰 度 ) FERRER 刀 年 轻 星 族 年 龄 
本 文 拟 合 值 0. 0250 6. 30 6. 30 
文献 [4] bsSP-fitted 拟 合 值 0.030+ 0.00 13.50 + 0.00 12.27 + 0.81 


分 析 产 生 误 差 的 原因 是 我 们 对 实验 模型 进行 了 简化 ， 一 是 理论 光谱 库 只 考虑 CZ, ty, t2) 而 
没有 考虑 v,, o TA4 等 参数 的 影响 ， 二 是 没有 考虑 观测 光谱 的 误差 .没有 对 观测 光谱 的 不 


确定 性 进行 处 理 ， 即 是 信 噶 比 (S/D 取 值 很 大 的 理想 情况 。 文 中 得 到 的 金属 丰 度 比 文献 [4] 


的 拟 合 值 小 0. 005。 这 也 符合 我 们 的 预期 ， 原 因 是 模型 输入 的 不 确定 性 通常 会 导致 金属 丰 度 


估计 值 偏 低 ”"。 星 族 


ER (ty, t2) 的 拟 合 值 比 文献 [4j 的 值 低 50% 的 主要 原因 是 简化 的 模型 


没有 考虑 速度 弥散 从 而 导致 星 族 年 龄 被 低估 。 该 拟 合 结果 也 体现 了 “年 龄 -金属 丰 度 简 并 效 


应 ” ` 


综 上 ,增加 自 


适应 策略 的 改进 遗传 算法 在 拟 合 速度 上 较 传 统 遗传 算法 和 BS2fit 算 法 平均 


提高 了 43. 5%， 拟 合 结 果 符 合 预 期 。 


4 总 结 与 讨论 


采用 自 适 应 策略 改进 的 遗传 算法 对 双星 星 族 光 谱 拟 合 ， 结 果 表 明 : d) 采用 群 智能 算 
法 (SIA) 解决 双星 光谱 拟 合 问题 是 可 行 的 。(2) 改进 策略 遗传 算法 和 传统 遗传 算法 、BS2fit 


算法 相 比 在 双 


星星 族 光 谱 拟 合 中 能 将 计算 速度 分 别提 升 19. 1% 和 67. 2%， 对 双星 星 族 光谱 拟 


合 速度 平均 可 提高 43. 5%， 提 升 效果 明显 。 改 进 策略 的 遗传 算法 在 传统 遗传 算法 的 基础 上 ， 


发 挥 遗 传 算法 全 局 搜索 能 力 强 、 相 较 于 传统 遗传 算法 收敛 速度 快 。 此 外 , 改进 策略 的 遗传 算 


法 在 搜索 时 更 


略 遗 传 算法 在 


智能 , 能 有 效 跳出 局 部 最 优 陷阱 ， 提 高 遗传 算法 全 局 搜索 能 力 , 避免 算法 早熟 ， 
相对 于 BS2fit 算法 的 全 网 格 搜索 SGA 算法 对 搜索 空间 做 了 取舍 所 以 效率 更 高 。 虽 然 改 进 策 


双星 星 族 拟 合 上 提升 了 速度 , 但 是 还 有 很 大 提升 空间 。 一 方面 ,编码 方式 对 算 


法 探索 能 力 和 收敛 性 也 有 影响 ,改进 算法 编码 如 利用 量子 编码 等 也 值得 探索 ， 另 一 方面 ， 群 


智能 算法 发 展 日 新 月 异 ， 探 索 更 有 效 的 搜索 策略 ， 如 探索 和 强化 学 习 Q-learning 等 融合 也 
很 有 必要 。 考虑 到 遗传 算法 易于 并 行 实现 , 下 一 步 可 探索 并 行 策略 下 的 改进 策略 遗传 算法 对 
光谱 拟 合 速度 的 改善 。 
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Binary-Star Spectral Fitting Based on Strategy 
Improved Genetic Algorithm 
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Abstract: The essence of the binary-star spectral fitting is to search the parameter space 


corresponding to the comparison of observed and theoretical spectra, via techniques such as 


least square. Rapid and accurate fitting of spectra is very important for studying a huge number 


of galaxies. Because it is not easy to deal with such problems as spectral fitting contains huge 
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amounts of calculation and the efficiency of fitting is not high. The strategic improved genetic 


algorithm is used for binary-star spectral fitting by this paper, and the comparison to BS2fit 


algorithm and traditional genetic algorithm is also shown. The experimental results show that 


the improved genetic algorithm can increase the spectral fitting speed by 43.5% on average, 


which can possibly promote the application of stellar population synthesis method in 


astronomical researches obviously. 


Key words: Spectral fitting; Galaxies; Improved Genetic Algorithm 


9/9 


